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1. Introduction 

The statistical properties of the spectra of dis- 
ordered systems continues to be a vivid field of re- 
search both analytically and numerically. It is by 
now well established that the spectral fluctuations 
in the limit of vanishing disorder are described by 
the random matrix theory (RMT) |l],|2j with cor- 
rections of the order of 1/g, where g is the dimen- 
sionless classical conductance of the system. This 
is because in this regime the states are extended 
(ergodic). In the limit of infinite disorder, on the 
other hand, due to the extremely strong localiza- 
tion of the eigenstates the levels are practically 
uncorrelated at least in the thermodynamic limit. 

The first attempts || to investigate the statis- 
tics at the localization-delocalization transition 
(LDT) (also known as the metal-insulator tran- 
sition), showed that indeed the statistics at the 
LDT is different from the two other extremes. 
Eventually it shares some of the features of both 
of them. It has been shown numerically in many 
simulations ||-|| that e.g. the nearest neighbor 
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spacing distribution, P(s), the probability of find- 
ing no eigenlevels within an energy window of 
length s, behaves for low-s as P{s) ~ s 13 , and for 
large-s as lnP(s) ~ —as, where f3 = 1,2 and 4 
depending on the global orthogonal, unitary and 
symplectic symmetry of the system, and a is a 
constant that depends weakly on j3 but strongly 
on dimensionality d The presence of the 

level repulsion at low-s is a direct consequence 
that the eigenstates at the LDT are heavily fluc- 
tuating but still strongly overlapping The 
heavy fluctuations expressed as their multifrac- 
tal character |l0[ enters in the parameter a that 
is connected both to the correlation dimension of 
the states and the level compressibility jflj . Such 
a relation shows how the statistical properties of 
the levels and the states are coupled at the LDT. 

In a previous publication Q we have presented 
numerical evidence that using appropriate indices 
to describe the overall shape of P{s), a supe- 
runiversal (i.e. [3 independent) scaling relation 
is revealed as the system undergoes the transi- 
tion from the metal to an insulator. Here we give 
further evidence that with the change of system 
size L and disorder W the nearest neighbor spac- 



ing distribution, -P(s), evolves as a one-parameter 
family of distributions. A similar evolution has 
been shown analytically to exist by Shapiro |T2|] 
for the conductance distribution in d = 2 + e. The 
parameter in our case is also the scaling variable 
that can be chosen as <?, i.e. the classical, dimen- 
sionless conductance of the system. 

2. The model and the analysis 

In order to study the statistical properties of 
disordered systems with and without time rever- 
sal symmetry and with broken rotational symme- 
try we have used the Anderson model which is a 
tight -binding model defined on a L 3 lattice. For 
the details we refer to the review of Kramer and 
MacKinnon []l3| and also our previous work ||. 
In Ref. (|] the numerical results have been illus- 
trated using the symplectic case where the rota- 
tional symmetry has been broken with the intro- 
duction of spin-orbit scattering. Here we present 
the results of the case with (without) time re- 
versal symmetry corresponding to the orthogonal 
(unitary) symmetry. Time reversal symmetry has 
been broken by the introduction of a magnetic 
field. 
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Figure 1. Q(L,W) and S(L, W) for the case 
of orthogonal symmetry. Continuous curves are 
polynomial fits. 
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Figure 2. Q(L, W) and S(L, W) for the case of 
unitary symmetry. Continuous curves are poly- 
nomial fits. 

After appropriately unfolding the spectra in 
many samples of the system we calculated the 
nearest neighbor spacing distribution, P{s). For 
each system size L and disorder W the number 
of spacings was of the order of 10 5 . In order to 
characterize the shape of P(s) we use generalized 
Renyi-entropics llJJ 
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where fi2 
while fis : 
obtain 

- ln(<?) 

Sstr 



(s 2 ) is the second moment of P(s), 
(sins). After a proper rescaling we 

- ln(q) + ln(q w ) ~ 



- bx(q F ) + ln(q w ) 

S s tr — S\Y 
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where index p refers to the Poissonian, uncorre- 
cted case and w to the Wigner-surmise repre- 
senting the RMT cases respectively. Their values 
are listed in Table [TJ. Such a rescaling ensures that 
§ = Q = 0(1) belonging to the RMT (Poisson) 
limit. In Figs. |l| and |^ we present the way how 
S and Q vary as a function of disorder for vari- 
ous system sizes. The scale invariance at around 
W — W c shows the presence of the LDT. 

In Table | we collected the values of the crit- 
ical disorder and the critical exponent obtained 
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Table 1 

Shape descriptive parameters for the case of dif- 
ferent P(s) functions. 

Poisson 0=1 = 2 = A 

q 0.5295 0.7854 0.8488 0.9054 

-ln(g) 0.6358 0.2416 0.1639 0.0994 

Sstr 0.2367 0.1025 0.0733 0.0464 

from a linearization around this fixed point. Even 
though the values presented in Table || agree well 
with the ones known from earlier calculations, the 
precision is somewhat lower. 

In Figure || we plot the S as a function of Q for 
all the symmetries ((3 = 1,2 and 4). As a com- 
parison we also plotted our analytical estimates 
based on a phenomenological assumption and the 
relation obtained for the interpolating formula of 
Israilev [p5f. We will show below evidence that 
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Figure 3. S(L, W) as a function of Q(L, W) 
for all the symmetry classes. The solid symbols 
represent the positions of the critical points. The 
continuous (solid, dashed, dotted) curves are our 
analytical estimates, see text for details. The 
dashed dotted line shows the relation of the P(s) 
due to Izrailev p5| representing the pure effect of 
level repulsion. 

the overall scaling behavior in Figure ^ can be un- 
derstood if the spacing distribution for any finite 



Table 2 

Position of the critical point and critical expo- 
nents obtained using the shape descriptive pa- 
rameters Q and S. 

W c v 

= 1(0) 16.77±0.63 1.30±0.38 

= 2(U) 18.13±0.38 1.63±0.57 

/3 = 4(S) 21.97±0.17 1.41±0.32 

value of W and L, i.e. for < g < oo is given as 

/oo 
Q 9 Ae X ~ V )me y )dy (3) 
-oo 

with e x = s. Here the function Wp(t) is simply 
the Wigner-surmise and V g ,p(s) is the histogram 
measured in the numerical simulation. All three 
functions, V, Q and W are spacing distribution 
functions. 

In the next section we will prove that for the 
convolution of the type (||) the Renyi-entropies 
introduced in are additive. Therefore, apart 
from the denominator, the constant shift in (Q) 
corresponds to Wp in (|J) , hence we may conclude 
that the relation in Figure || is described by the 
properties of function Q. 

We will also show that the solution of (||) for 
Q,g t fj(u) in the extreme cases (g = and g — > oo) 
can be done analytically and that a phenomeno- 
logical interpolation can describe the S vs. Q 
relation in Figure ||. 

3. Additivity of Renyi entropies 

The expression ([}]) is in fact a convolution and 
after a variable transformation can be written as 

/"DO 7i 

v g A*) = J o Q B ,p(t)Wfi®j' (4) 

Let us denote the k-th moment of a distribution T 
by Hk [F] ■ Then it is easy to see from (||) (omitting 
indices g and 0) that 

/>oo 

Hk[V] = / dss k V(s) 

- fe^r***® 

poo poo 

= / dtt k W{t) I duu k Q(u) 
Jo Jo 
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= MWW[Q], (5) 

which for q = — In /i2 yields the desired additivity. 
A similar procedure yields the additivity of the 
quantity S s t r - The key step is to realize that the 
entropic moment of T is ^s[.F] — (sins) which 
can be expressed as 

ns[P] = MQ]ms[W] + m\W\ns[Q]- (6) 

4. Phenomenological approximation 

Here we show that the solution of (||) for 
Qg,p(u) in the extreme cases (g = and g — > oo) 
can be obtained analytically and that a phe- 
nomenological interpolation can describe the 5 
vs. Q relation in Figure |[ First it is easy to see 
that for the case of vanishing disorder (g — > oo) 
we expect that V{s) — ► W(s). In this case 
<2(u) = <5(u — 1). On the other hand at the other 
extreme, when g — > 0, another universality is ex- 
pected: V{s) — > exp(— s), the Poissonian distri- 
bution. In this case after nontrivial calculations, 
Eq. (g) is solved for Q ,p(s) = R{bs) 

f e-y 2 (O) 
fl(y)=<J erfc(y) (U) (7) 

[ (2y 2 + l)erfc(y) - (S) 

where y — bs, with b = 1/t/tt, \/tt/4, and 3y/n/16 
and a = 2/tt, tt/4, and 9tt/32 for /3 = 1, 2 and 4, 
respectively. The complementary error function 
is denoted as erfc. Now a phenomenological step 
is introduced in order to describe the evolution of 
V for any finite value of g. We propose a simple 
form 

Q g As) = a g ,^R(b g<0 s) (8) 

where the coefficients a ffl/ g and b g> p are deter- 
mined from the normalization conditions /io[Q] = 
Mi [2] = 1- Obviously for g — > they yield 
the values given above for (Q). For g — > oo 
<2(s) — > <5(s — 1), as well. The 5 vs. Q relation 
for this interpolating function is depicted as con- 
tinuous curves in Fig. ||. They still have a small j3 
dependence but follow the data very closely. The 
/3-dcpcndence can be attributed to (i) the simple 
interpolation we used in (||) and to (ii) the fact 
that we focused only on the linear shift introduced 
in (Q) and disregarded the denominator. 



5. Conclusion 

In this paper we have given further details, 
both numerical and analytical of a recently dis- 
covered scaling relation existing in the evolu- 
tion of the nearest neighbor spacing distribution, 
P(s), as the dimensionless conductance g changes 
from (localized) to oo (metal). 
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